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Abstract In finance, the price of a volatile asset can be modeled using frac- 
tional Brownian motion (ffim) with Hurst parameter H > \. The Black- 
Scholes model for the values of returns of an asset using fBm is given as, 

Y t = Y exp ((r + n)t + oB? ) , t > 

where Yq is the initial value, r is constant interest rate, y, is constant drift and 
a is constant diffusion coefficient of fBm, which is denoted by (B^) where 
t > 0. Black-Scholes model can be constructed with some Markov processes 
such as Brownian motion. The advantage of modeling with fBm to Markov 
proccesses is its capability of exposing the dependence between returns. The 
real life data for a volatile asset display long-range dependence property. For 
this reason, using fBm is a more realistic model compared to Markov processes. 
Investors would be interested in any kind of information on the risk in order 
to manage it or hedge it. The maximum possible loss is one way to measure 
highest possible risk. Therefore, it is an important variable for investors. In 
our study, we give some theoretical bounds on the distribution of maximum 
possible loss of fBm. We provide both asymptotical and strong estimates for 
the tail probability of maximum loss of standard fBm and fBm with drift and 
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diffusion coefficients. In the investment point of view, these results explain, 
how large values of possible loss behave and its bounds. 

Keywords First keyword • Second keyword • More 
1 Introduction 

In finance, the price of one share of the risky asset, is modeled using fractional 
Brownian motion (fBm) with Hurst parameter H € [1/2, 1) in order to display 
long-range dependence. Hence, our results are on fractional Brownian motion 
with Hurst parameter H € [1/2, 1). 

There are several stochastic integral representations that have been devel- 
oped for fBm. For example, 

»3w^i ((t -' )H "^- ( -' )H " 1P)d,,r - 

+ / (t~s) H -^ 2 dW s ) (1) 
Jo 

is a fBm with Hurst parameter H £ (0, 1), where W is a Wiener process. 

Let H be a constant in the interval (0, 1). A (standard) fBm {Bf : t > 0} 
with Hurst parameter H is a continuous and centered Gaussian process with 
covariance function 

E[B»B?] = \{t™ + s™-\t-sr) (2) 

and Bq — 0. It follows that B H has stationary increments, that is B^ +s — Bf 
has the same law as B^ , for s,t > 0. 

For H = 1/2, the process (B^)t>o corresponds to a standard Brownian 
motion, in which the increments are independent. By definition, the covariance 
between the increments B H (t+h)-B H (t) and B H (s+h)-B H (t) with s+h < t 
and t — s — nh is 

p H {n) = \h™[{n + lf H + {n~ lf H - 2n™]. 

We observe that two increments of the form B H (t + h) — B H (t) and B H (t + 
2h) — B H (t + h) are positively correlated for H > 1/2, and they are negatively 
correlated for H < 1/2. 

Since the covariance function of fBm is homogeneous of order 2H, [5] fBm 
possesses the self-similarity property, that is, for any constant c > 0, 

{B ct ) t >o - (C B t )t> a . 

The aim of the present paper is to find the distribution of maximum possi- 
ble loss of fBm, which is one way of measuring investors risk. In order to obtain 
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results related to this distribution, we are also interested in the supremum, the 
infemum, and the range variables of fBm. 

Let (i?, 9, P) be a probability space and let B H be fractional Brownian 
motion with Hurst parameter He [1/2, f). We introduce the following nota- 
tion. 

— Let I? := inf B„ denote the infemum of fractional Brownian motion up 

* 0<v<t v 

to time t. 

— Let := sup B v denote the supremum of fractional Brownian motion 

0<v<t 

up to time t. 

— Let — ~ , called the range of fractional Brownian motion up to 
time t. 

— The maximum loss of fractional Brownian motion before time t is defined 

as 

M? := sup (B u - B v ) = sup ( sup {B u - B v )) 

0<u<v<t 0<v<t 0<u<v 

More generally, the above definitions can be repeated for a non-standard 
fBm. Let \x be drift parameter taking real values other than and similarly 
a > be a real valued diffusion coefficient other than 1 for fBm with drift 
defined as Y t :=fit + aB t . 



2 Bounds on the distribution of maximum loss 

In this section we introduce some new bounds on the expected value of max- 
imum loss of fBm and on the distribution of maximum loss of fBm. These 
results are already useful for investors as they are. Also, in later sections they 
will be useful for finding the asymptotic distribution of maximum loss. 

Theorem 1 For fBm up to time a with Hurst parameter H > i, and for 
y > 0, we have 

2V 7r V 71 " 
and _ 

P(Mf - > y) < P{R* > y) < 2 -^f- 

Proof: E(S^) < ^2 * a H is known 3 and by the Markov's inequality an 
upper bound for the distribution of the supremum is found, that is for x > 0, 
P(S^ > x) < ■ Then, by the symmetry property of centered Gaussian 

processes one can show that E(I^) > — ^= * a H 

Combining the results given above we find an upper bound for the expected 
value of range, Rf, that is E(R^) < '^j=- * a H . Furthermore, by Markov's 

inequality we see that for y > 0, P(Rf > y) < 
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Clearly one can see that, 

Jf := - inf X v < sup sup (X ( f - X?) = M t H - < R? 

0<v<t u<v<tO<u<v 

holds for all t. 
Hence 

- E{I?) < E(M?>-) < E(R») < 2 -^f- (3) 
is obtained, and by Markov's inequality we get 

P(Mf - > y) < P(R» >y)< ^f- 

We have also noticed that, < E(S?) < These bounds 

are obtained using H = 1 and H = 1/2 in Sudakov-Fernique inequality [U 
Theorem II.2.9]. 

Furthermore, because E(S^) equals —E(I^) using Equation ([3]) we obtain, 

^ < E(M «n < 

2V7T 

□ 



3 Asymptotical distribution of maximum loss of fractional 
Brownian motion 

The theorem in this section gives the asymptotic result for the probability 
distribution of maximum loss of fBm. We first start with scaling property of 
loss process which we denote as X v := sup (B u — B v ), v > 0. Our approach 

0<u<v 

is similar to the large deviations technique used for queeing systems modeled 
by fractional Brownian motion, [5] and [5]. 

We define the loss process X by X v = sup (B u — B v ). 

0<u<v 

Proposition 1 The loss process X is self-similar and each X v has the same 
distribution as the supremum S v for every v > 0. 

Proof: The self-similarity of fBm corresponds to {B au : u > 0} ={a H B u : u > 
0} for every a > 0. It follows that 

{B au - B av : < u < v,v > 0} = {a H (B u - B v ) : < u < v, v > 0} . 
Therefore, we get 

{X av : v > 0} = {a H X v : v > 0} 
by the definition X v = sup (B u — B v ). 

0<u<v 
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On the other hand, since fractional Brownian motion has stationary incre- 
ments, the collections {B u — B v : < u < v} and {— B v - U : < u < v} have 
the same probability law for fixed v. Both are mean Gaussian processes with 
covariance function 



r(u, u') = 1/2 [ \v - u\ 2H + \v- u'\ 2H - |- - n2H i 



— \u — u i 

Since the suprcmum of the two collections will also have the same distribution 
and {-B v _ u : < u < v} ={B U :0<u<v},we get X v = . □ 

Now, let us use # to denote cumulative distribution function of Normal 
distribution and <P for its complement distribution 

Proposition 2 For all x € R + , we have 

P{M t H '~ >x)> ${x/t H ). 

Proof: As a result of Proposition [TJ for < v < t, we see that, 

P{X V >x)=P( sup B u > x) > sup P{B U > x) = sup 3{x/u H ) = ${x/v l 

0<u<v 0<u<v 0<u<v 

(4) 

Similarly, we observe 

P(Mf'~ > x) = P( sup X v > x)> sup P(X V >x)= sup 3(x/v H ) = $(x/. 
o<v<t a<v<t o<v<t 

(5) 

Theorem 2 For the maximum loss M t ' of fBm and x > 0, we have 

Proof: We start with finding a lower bound for limit internum of the logarithm 
of the distribution of maximum loss. Combining Proposition [5] with 

lim {l/x 2 )\og${x/t H ) = -{2t 2H y 1 
(Adler, 1990, sf.42), we obtain, 

Uminf \ \ogP{Nl"'~ > x) > lim \ \og$(xlt H ) = --^7 (6) 

x— >oo x x— >oo X 2t 

We continue with finding an upper bound for the limit supremum. Note that 
maximum loss M^'~ has two dimensional index set T :— {(u, v) : < u < 
v < t}, and it is a centered Gaussian process. Based on the fact that T is a 
separable metric space we use Borel's inequality given in jX, Theorem II. 2.1] 
for finding the supremum of this process. By the continuity property of the 
paths of fractional Brownian motion, B U ~B V , is bounded on T. And by Borel's 

inequality or directly from Theorem[T]we see that r\ :— E [ sup (B u — B v )\ 

a<u<v<t 

is finite. 
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As a result of Borel's inequality 

P(M t H >x) = P{ sup (B u - B v ) >x)< 2e-^ x -' n)2 /t2 " , x > rj 

0<u<v<t 

can be written. Here, sup E(B U — B v ) 2 — t 2H is used. 

0<u<v<t 

This shows, 

lim sup — log P(M^'~ > x) < lira 

x — s-oo X x—>oo 

which completes the proof. □ 

In [B] it was shown that the result given in Theorem [5] for centered, Gaus- 
sian sequence of random variables which was mentioned in [TJ pg.43]. How- 
ever, in order for us to use this result, it must be shown for general metric 
space T, similar to the proof of Borel's inequality. In stead of showing this 
generalization, we preferred giving the proof for M t '~. As a conclusion, we 
see that the increase in the logarithm of distribution of M t '~ behaves same 
as the asymptotic increase of distribution of B t . Here, the important property 
of B t is, it has the same distribution as the random variable with the highest 
variation among the random variables {B u — B v : < u < v < t} which is 
Bq — B t = —Bt- 



log 2 1 (x — rjY 



2t 2H 



(7) 



4 Asymptotical Estimate for the Tail Probability 

In the next theorem, we generalize the previous result for fractional Brownian 
motion Y with drift and diffusion coefficient. Recall Y t = [it + aB t for /i 6 
R, a > 0. 

Theorem 3 The maximum loss process defined by M^' := sup (Y u — Y v ) 

0<u<v<t 

satisfies 

1 , , 1 



}™^2 ^ P (M t > >x) = - 2aH2H 
where Y is a fractional Brownian motion with drift and diffusion coefficients. 
Proof: For all x E the trivial lower bound is given by 

P(M t H <~ > x) = P( sup (Y u -Y v )>x)> sup P(Y U -Y v >x). 

0<u<v<t 0<u<v<t 

Since Y has stationary increments, we have 

P(Y U -Y v >x) = P(-Y v _ u >x) = 3((x + fi(v - u))/(a(v - u) H )) . 
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Therefore, the following hold 

„, . = / X + U(v — u)\ - ( X + LLv\ 

sup P(Y U -Y v >x)= sup <M p J = sup # £- 

o<m<«<* o<w<ti<t V a(v — u) J o<v<t \ crv j 

Since ^ is a decreasing function, we find v G [0, t] that minimizes f(v) = 
(x + fiv)/(av H ). The critical value of / is v* — j^j^m ■ We check if v* € [0, t] 
or not as below 

— If /i > and x < t/i(l — H)/H, then the minimum value of / is obtained 
at v* £ [0,t], 

— If > and x > tfi(l — H)/H, then the minumum of / is obtained at t. 

— If /j. < 0, then the minimum of / occurs at t. 

As a result, we have P(M^'~ > x) > $((x + fit)/(at H j) for large x > 0. 
We get 

liminf \ log P(Mf<- > x) > Urn L\ g§{^±£) = -— ^ . 

On the other hand, we use Borel inequality which applies to a mean zero 
Gaussian process in order to find the limit supremum [I] Theorem II. 2.1]. By 
definition, Y u — Y v = a(B u — B v ) + pb(u — v) and 

sup (Y u - Y v ) > x 3(u, v) G t , a(B u - B v ) + fi(u - v) > x 

0<u<v<t 

cr(B u -B v ) 
3(u, v)et, — V -A- > l 

x — fi(u — v) 
o(B u - B v ) 

sup > 1 

0<u<v<t X - fi(U - V) 

where we assume x > [it. Making a change of variable (u, v) — > (u/x, v/x) and 
using the self-similarity of standard fractional Brownian motion, we get 

P( sup (Y u -Y v )>x)=p( sup < B ^-B V ) 

0<u<v<t \0<u<v<t X — /J,(U — V ) 

tj I cr(B xu - B xv ) \ 
= P sup > 1 

y0<«<1;<| X-fIX{U-V) J 

( ax H (B u -B v ) \ 
= P\ sup — r > x \ . 

Since G u v :— a (B u — B v ) / (1 — jj,(u — v)) is a zero-mean continuous Gaussian 
process, Borel inequality implies that for x > rj\ 



P(M t H '~ > x) = P{ sup (Y„-Yj>x)<2exp 

0<u<v<t 



1 {x 1 '" - r) x ) 



, 2 "1 



(8) 
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where rj x = E [ sup G u v ], 7 = sup E G 2 U v , and it is observed that 

0<ti<u<t/a: 0<u<v<t/x 

Vx < ??i holds for large x, in particular x > t. We have 

sup & = sup — = sup — 

0<u<v<t/x 0<u<v<t/x 1 + ~ u ) 0<v<t/x l +l lv 



and therefore 



< Tj x < sup 



E{ sup By) 



(9) 



0<v<t/x 1 + ^ V 0<v<t/x 

Note that (t/(1 + jiv) is bounded from above by a for large x. Since 

V^(t/x) H V2(t/xf 
— -= — < E{ sup B v ) < r= , 

0<t><t/z v 71 " 

we take i] x oc (t/x) H below in view of ©. On the other hand, when x is large 
enough, that is, when tjx < H/(l — /iH), 

° 2 v 2H a 2 (t/x) 2H 
7 = sup EG UV = sup — -2 = 2 

0<u<tJ<t/ic 0<v<t/x [1 + fiV) [1+fXt/X) 

is found. As a result of the bound ©, we get 



limsup — log P(M^' > x) < lim 



in view of H G (0, 1) 



1 (a; 1 -" - 



2a 2 {t/x) 



2H 



1 



□ 



5 Stronger Form of the Asymptotic Distribution 

In this section, we directly find the asymptotical form of the tail probability 
using a characterization of [3] for Gaussian processes. 

Theorem 4 Asymptotically, we have 

,. JP(M t H - > x) 

hm — , , = 1 . 

x^oo <P(x/t 2H ) 

when M t is defined for standard fractional Brownian motion. 

Proof: The proof is based on two conditions given in [5] that characterize the 
existence of the limiting distribution. We work with the mean zero Gaussian 
process {B u — B v : < u < v < t}. Let a\ = sup( u t ,) gT JE(i?„ — B v ) 2 , which 
yields Oy, = t 2H . The first condition is that there exists a unique (uq,vq) in 
T such that JE(B Uo - B Vo ) 2 = a\. This holds with (u ,v ) = (0,t). For the 
second condition, the set 



T h = {{u,v) 6 T : M(B t (B u - B v )) > a 2 T - h 2 } 
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is defined for ft, > 0. In order to identify Th, we find W l (B t (B u — B v )) — 
\{v 2H - u 2H + (t - u) 2H -(t- v) 2H ). Therefore, (u, v) G T H satisfy 

y 2H _ U 2H + (t _ u] 2H _ {t _ v) 2H > rfH _ 2/j 2 (1Q) 

This implies that (fTU|) is satisfied by (u, v) G Tff, for fixed u. Now, for u > t/2, 
we have (t - u) 2H - u 2H < 0, and (t - v) 2H > for all v G [0, <]. Then, from 

m 

y 2H > 2i 2H _ 2/l 2 > j2H _ 2 ^2 

On the other hand, for u < t/2, we have it — u) 2H < t 2H . Therefore, we get 

v 2H + t 2H > v 2H - u 2H + (t - u) 2H - (t - t;) 2H > 2< 2ff - 2ft 2 

which again implies t> 2ff > i 2ff — 2ft 2 for fixed u and (fi, v) G T/j. Since, 
/(d) = i; 2ff is convex, we get 

v>t- Kh 2 

for some constant K [5J pg.309]. Now, we consider second condition in [5] 
which requires 

lim ft" 1 IE sup B u -B v +B t = l (11) 

In particular, we have 

IE sup B a - B v + B t < IE sup B a - B v + B t 

(u,v)eT h v>u, t-v<Kh 2 

Then, it follows that 

lim ft" 1 IE sup Bu- B v + B t = lim ft" 1 IE sup B t - V (12) 

v>u : t-v<Kh 2 t-v<Kh 2 

since fractional Brownian motion has stationary increments and TEB a = 0. 
On the right hand side of (|12p , the supremum of fractional Brownian motion 
[0, Kh 2 } is bounded by y/2K H h 2H /y^F [3[5] and hence we get the limit in (12]) 
to be 0. Since T is separable, a monotone convergence argument extends the 
result for fixed u to all Th proving (fTTj) [TJ pg.47]. □ 
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